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Research  Results 


This  description  of  results  from  the  project  covers  research  in  two  areas: 

1.  Plastic  flow  in  two  and  three  dimensions. 

2.  Hyperbolic  conservation  laws. 


1.  Plastic  flow.  This  part  of  the  project  was  carried  out  in  collaboration 
with  David  Schaeffer  of  Duke  University.  Our  research  focused  on  the  issue 
of  loss  of  stability  and  well  posedness  in  the  equations  of  motion  of  granular 
materials.  The  partial  differential  equations  are  derived  from  conservation 
of  mass  and  momentum,  augmented  by  constitutive  laws  that  relate  the 
dependent  variables  algebraically.  The  starting  point  was  the  motion  in  two 
dimensions  of  a  rigid-plastic  material,  with  the  constitutive  laws  coming  from 
critical  state  soil  mechanics.  We  established  a  criterion  for  stability  in  [8],  and 
related  this  to  the  quasidynamic  equations  in  [9],  showing  broadly  speaking 
that  the  full  equations  are  linearly  stable  if  and  only  if  the  quasidynamic 
equations  are.  This  is  important  because  the  quasidynamic  equations  are 
much  more  convenient  for  computation,  and  our  results  show  that,  with 
precisely  formulated  conditions,  the  simpler  equations  capture  the  stability 
properties  of  the  original  system.  A  third  paper  in  preparation  analyzes 
corresponding  issues  for  flow  in  three  dimensions. 

We  recently  looked  at  yield  vertex  constitutive  equations.  These  differ 
markedly  from  those  of  critical  state  soil  mechanics,  in  that  the  equations 
become  fully  nonlinear  rather  than  merely  quasilinear.  We  found  that  the 
equations  lose  hyperbolicity,  a  phenomenon  mathematically  associated  with 
a  loss  of  well-posedness,  and  physically  associated  with  localization,  leading 
to  the  failure  of  the  material  through  the  formation  of  shear  bands.  A  no¬ 
table  success  of  the  analysis  in  [10]  is  the  correct  prediction  of  the  preferred 
orientation  of  the  shear  band  in  an  experiment  in  which  a  granular  material 
is  sheared  between  parallel  plates. 

2.  Hyperbolic  conservation  laws.  The  classification  [A]  of  2x2  sys¬ 
tems  of  hyperbolic  conservation  laws  with  quadratic  nonlinearities  identifies 
four  different  types  of  equations.  The  Riemann  problem  was  solved  in  detail 


in  [B]  for  three  of  the  four  types.  The  fourth  type  of  equation,  Case  I,  is  the 
most  significant  for  applications  to  models  of  multiphase  flow  in  oil  reser¬ 
voirs,  as  discussed  in  [A].  This  case  involves  undercompressive  shocks,  which 
are  physical  shock  waves  closely  associated  with  systems  that  change  type. 

Using  ideas  from  dynamical  systems  to  understand  the  role  of  undercom¬ 
pressive  shocks,  the  Riemann  problem  was  solved  for  Case  I  equations  in 
[1].  This  paper  exploits  the  special  properties  of  quadratic  nonlinearities.  In 
particular  the  solution  of  the  Riemann  problem  is  not  stable  to  perturbation 
of  the  nonlinearity  by  higher  order  terms.  This  problem  was  solved  using 
equilibrium  bifurcation  theory  augmented  by  a  characterization  of  hetero- 
clinic  orbits  using  Melnikov’s  integral.  Preliminary  results  were  announced 
in  conference  proceedings  [2,3],  and  written  up  in  detail  in  [4]. 

A  detailed  study  of  change  of  type  for  model  equations  of  three  phase 
flow  in  porous  media  is  given  in  [6].  These  equations  typically  have  small 
elliptic  regions  that  may  be  somewhat  accentuated  by  the  inclusion  of  gravity 
effects,  as  noted  in  numerical  results  with  specific  models.  Special  classes  of 
equation  are  now  known  to  lose  strict  hyperbolicity  while  remaining  hyper¬ 
bolic.  Corners  of  the  physical  domain  correspond  to  umbilic  points,  at  which 
the  characteristic  speeds  coincide,  but  they  are  degenerate,  leading  us  to  an 
analysis  of  the  higher  order  terms. 

There  are  classes  of  nonstrictly  hyperbolic  systems  for  which  the  char¬ 
acteristic  speeds  are  real,  but  coincide  along  a  curve.  In  [7],  we  give  a  lo¬ 
cal  analysis  and  classification  of  such  equations,  and  show  how  new  types 
of  shocks,  known  as  singular  shocks,  are  a  crucial  part  of  solving  Riemann 
problems.  We  give  a  detailed  asymptotic  analysis  of  corresponding  solutions 
of  the  regularized,  parabolic,  equations.  These  solutions  blow  up  at  a  single 
point  as  the  dissipation  approaches  zero. 

Glimm’s  method  was  implemented  for  the  full  initial  boundary  value  prob¬ 
lem  describing  the  motion  of  an  elastic  string  stretched  between  two  fixed 
points  [1 1] .  The  numerical  results  are  strongly  indicative  of  the  presence  of  a 
periodic  solution.  This  is  reinforced  by  the  discovery  of  two  exact  solutions 
that  describe  the  motion  of  sections  of  the  string.  The  analytic  solutions 
can  be  combined  to  form  a  periodic  function  that  resembles  the  numerical 
solution. 
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